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METHOD AND APPARATUS FOR POLARIZATION 
CONTROL WITH PHOTONIC CRYSTALS 



Cross-Reference to Related Applications 

This application claims the benefit of and priority from United States Provisional 
Application Serial No. 60/442,999 filed on January 29, 2003 and United States Provisional 
Application Serial No. 60/505,046 filed on September 24, 2003. 
Introduction. Background and Summary of Invention 

The present invention relates to polarization control in general. More particularly, the 
invention provides for the first time, a method and apparatus for complete polarization 
control of electromagnetic (EM) waves using photonic crystals. The novel apparatus 
utilizing the invention includes, for example, photonic-crystal-based transmission 
waveplates, polarizers, and reflection waveplates. These devices have tremendous 
advantages over the prior art: they can be much smaller than prior art devices, of the order 
of magnitude of one thousand times smaller; they can also be adjusted and modified to fit 
end uses more extensively than the prior art; and they can be scaled to any wavelength 
range including, but not limited, to the microwave, infrared, visible, and ultraviolet portions 
of the EM spectrum. It is also important that the birefringent properties of photonic crystals 
are related to the lattice geometry and indices of refraction of the constituent materials; 
thus, the birefringent and/or dispersive properties of photonic crystals can be tailored by 
adjusting the lattice geometry and size, as well as altering the constituent materials and 
their relative proportions. 

Photonic crystal birefringence has a range of potential applications. For example, these 
and related photonic-crystal devices will be useful in laboratory equipment, and to control 
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the flow of light in circuits and all-optical devices for communication and computing 
purposes. 

Furthermore, the use of photonic crystal birefringence and dispersion can be utilized in 
nonlinear optics to achieve phase-matching. This will afford many advantages over 
conventional phase-matching techniques. For example, in the prior art, phase matching in 
nonlinear optical processes is typically accomplished with birefringent, nonlinear materials. 
Although there are many types of materials with optical nonlinearity, few are also 
sufficiently birefringent. However, a photonic crystal can be birefringent and/or dispersive 
even if its constituent materials are neither birefringent, nor dispersive; therefore, using a 
photonic crystal constructed from a material with optical nonlinearity, it is possible to 
achieve phase matching without a material that is intrinsically birefringent. 

It is also possible to utilize photonic-crystal birefringence to couple two different modes 
(e.g., modes of polarization) for the purpose of creating an optical delay line. This type of 
delay line can in fact be created using any birefringent medium, which need not 
necessarily be photonic-crystal-based. 

A primary object of the invention is to provide a method and apparatus for controlling 
and/or manipulating the polarization of electromagnetic (EM) waves with photonic crystals. 

A further object is to provide a photonic crystal transmission waveplate capable of 
manipulating the polarization of light transmitted through the photonic crystal. 

A further object is to provide a photonic crystal reflection waveplate capable of 
manipulating the polarization of light reflected by the crystal. 

Another object of the invention is to provide a photonic crystal polarizer capable of 
reflecting one polarization and transmitting the other polarization (of a single wavelength), 
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and which is not limited by Brewster's angle, so that the polarized reflected and 
transmitted beams can have any angle relative to the incident beam. 

A further object is to provide a method of polarization control using photonic crystals for 
controlling the flow of light in circuits and all-optical devices used in laboratory equipment, 
communications, and computers. 

Another object is to provide a method of polarization control using photonic crystals to 
achieve phase-matching in nonlinear optics. 

Another object is to provide a method of creating a novel delay line for electromagnetic 
(EM) waves using energy transfer between different electromagnetic modes. 

Other objects and advantages of the invention will become apparent from the following 
description and drawings wherein: 
Brief Description of the Drawings 

Fig. 1A is a top or plan view of an example of a two-dimensional (2D) photonic crystal 
used in the present invention wherein the periodic plane or plane of incidence is parallel 
with the plane of the drawing; 

Fig. 1B is a side perspective view of the two-dimensional (2D) photonic crystal 
illustrated in Fig. 1A and wherein the optic axis is shown as A-A and extends in a direction 
parallel to the longitudinal axis of each of the elongated (extruded) hollow rods which form 
the illustrated crystal; 

Fig. 2A illustrates a photonic crystal quarter-waveplate according to the present 
invention and shows an incident linearly polarized beam entering the crystal and being 
transformed to a circularly polarized wave by the crystal; 
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Fig. 2B is a schematic representation of a photonic crystal half-waveplate utilizing the 
present invention wherein orthogonal components of an incoming polarized wave 
experience a relative phase shift of 180° as they are transmitted through the crystal; thus, 
the polarization of an incoming beam is rotated by the crystal as shown; 

Fig. 2C is a schematic representation of a photonic crystal polarizer of the present 
invention utilized to transmit one specific polarization through the crystal and 
simultaneously reflect the other polarization; 

Fig. 3 is a schematic representation of a photonic crystal slab showing the slab in a 
side sectional view; 

Figs. 4A and 4B are graphical representations of experimentally measured indices of 
refraction in the transparent spectral regions of hexagonal photonic crystals (lattice 
constant 1/2 inch, air-filling fraction [AFF] 0.60); 

Figs. 4C and 4D are graphical representations of the experimentally measured 
differences between the TE and TM indices of refraction in the transparent spectral 
regions of hexagonal-lattice photonic crystals, along with corresponding half-wave 
conditions (HWC) and quarter-wave conditions (QWC); 

Fig. 5A is a graphical representation of experimental measurements of intensity 
transmission through a photonic crystal as a function of frequency; 

Fig. 5B is a graphical representation of experimentally measured, normalized field 
(amplitude) transmission versus rotation angle for a hexagonal-lattice, photonic crystal 
quarter-waveplate (horns aligned - filled squares, horns crossed - open circles); 
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Fig. 5C is a graphical representation of experimentally measured, normalized intensity 
(power) transmission versus rotation angle for a hexagonal-lattice, photonic crystal half- 
waveplate (horns aligned - open circles, horns crossed - filled squares); 

Figs. 6A and 6B are graphical representations of experimentally measured normalized 
field (amplitude) transmission through a 20-layer, hexagonal-lattice photonic crystal over 
specified frequency ranges (horns aligned - open circles, horns crossed - filled squares); 

Fig. 7 is a graphical representation illustrating the existence of birefringence in a 
photonic crystal with respect to frequency (based on a simple model); 

Fig. 8 is a schematic representation of an experimental setup for conducting a variety 
of photonic crystal experiments; 

Fig. 9 is a graphical representation of normalized intensity (power) reflectivity in the 
vicinity of the fundamental band gaps of a 30-layer, square-lattice photonic crystal; 

Fig. 10 is a graphical representation of the experimentally measured phase of waves 
reflected from a 30-layer, bulk, square-lattice photonic crystal (lattice constant 1/2 inch, air- 
filling fraction [AFF] 0.66) in the vicinity of the fundamental band gaps as a function of 
frequency; 

Fig. 1 1 is a graphical representation of the experimentally measured reflected field 
(amplitude) versus receiver horn angle (open circles represent horns aligned, and solid 
squares represent horns crossed); these results were achieved using a 30-layer, bulk, 
square-lattice photonic crystal at approximately 9.56 GHz, ideal quarter-waveplate curves 
are shown in dotted (horns aligned) and solid (horns crossed) lines; 

Fig. 12 is a graphical representation of the experimentally measured difference 
between the TE and TM indices of refraction for a 20-layer, bulk, hexagonal-lattice 
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photonic crystal (lattice constant 1/4 inch, air-filling fraction [AFF] 0.32) well below the 
fundamental band gaps; also shown as a dotted line is the (f order quarter waveplate 
condition; 

Fig. 13 is a graphical representation of the calculated effective indices of refraction for 
10-, 20-, 40-, 60-, 80-, and 100-layer, one-dimensional (1D), acrylic-air, dielectric stacks of 
the indicated air-filling fractions (AFFs); the 0.62 AFF stack is a quarter-wave stack while 
the other two were created by both reducing the air layer thickness and increasing the 
acrylic layer thickness, or vice versa; the calculations were performed with the transfer 
matrix method for dielectric slabs; 

Fig. 14 is a schematic representation of a nonlinear photonic crystal used for second 
harmonic generation; 

Fig. 15 is a schematic representation of an experimental setup for demonstrating a 
photonic-crystal-based, mode-interference delay line; and 

Figs. 16A and 16B are graphic representations showing calculated (dotted lines) and 
experimentally measured (solid lines) transmission on the vertical axis and phase shown 
on the horizontal axis for the detected polarization inclined at 40° in Fig. 16A and inclined 
at 50° in Fig. 16B. 

Detailed Description 

Theory of Operation and Experimental Results 

A photonic crystal as known in the prior art is a periodic crystalline structure that 
possesses a unique band structure for electromagnetic (EM) waves. Typically, photonic 
crystals are composed of dielectric materials arranged in an array with a well-defined 
geometry, periodic in one or more dimensions. EM waves incident on the structure 
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produce multiple Bragg reflections from the many material interfaces that compose the 
structure. These multiple reflections interfere with each other to produce band gaps in the 
spectrum of the structure. Radiation whose frequency and wavevector lies within a band 
gap is exponentially attenuated as it propagates through the crystal. The characteristic 
band structure of a photonic crystal is specific to its geometry and material composition, 
including the indices of refraction and relative proportions of the different components. 
Although the basic physics behind photonic structures is application of Maxwell's 
equations, the complex boundary conditions at the material interfaces produce a wealth of 
interesting and novel phenomena. Since Maxwell's equations are scale invariant, results 
obtained at any wavelength apply to photonic crystals designed for any portion of the EM 
spectrum. 

There are many useful analogies between the EM band structure of photonic crystals 
and the electronic band structure of solid-state materials. However, there is an essential 
difference between the two cases: scalar theory is usually adequate in solid-state physics, 
but a full vectorial approach is in general necessary for photonic crystal calculations. This 
difference has been largely overlooked, leaving many interesting polarization-dependent 
phenomena of the present invention heretofore undiscovered. 

Any photonic structure with anisotropy should in general show polarization-dependent 
properties. In particular, if different polarizations experience different boundary conditions, 
the crystal should display a polarization-dependent band structure. 

In addition to the well-known power transmission effects of photonic crystals, nontrivial 
phase delays are also imparted to waves traveling within them. Although the power or 
intensity effects of photonic crystals have been well studied, phase effects, especially in 
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transparent spectral regions, have been heretofore ignored. It is possible to understand 
these phase effects by analogy with absorptive systems. Although dielectric photonic 
crystals have negligible absorptive dissipation, they have frequency-dependent reflectivity. 
In particular, energy that is not transmitted by the band gap is instead reflected. It is well 
known that a spectral region of absorption generates effects on the phase of transmitted 
waves, even in spectral regions far from the absorption line (i.e., transparent regions). 
This phenomenon is in fact required by causality through the Kramers-Kronig relations. 
Since the Kramers-Kronig relations are not specific to absorption only, we can expect 
similar results if they are applied to the frequency-dependent transmission profile of a 
photonic crystal band-structure. Thus, if the band structure of a photonic crystal is 
polarization-dependent, application of the Kramers-Kronig relations shows that there must 
be polarization-dependent (i.e., birefringent) phase effects, even in transparent spectral 
regions. 

The polarization-dependent properties of photonic crystals are, for the first time, 
exploited with the present invention for the purpose of controlling the polarization of EM 
waves. In particular, photonic crystals are used for the first time as novel, compact 
waveplates, polarizers, and polarizing-beamsplitters capable of polarization control that 
can be precisely tailored by modifying the material composition as well as the geometry of 
the structure. 

Photonic Crystals Used in Experiments 

As practical examples of these devices, it is useful to consider and conduct 
experiments with two-dimensional (2D) photonic structures as shown in Figs. 1A and 1B. 
These structures have a 2D-periodic lattice and are non-periodic or extruded in their third 
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dimension. Bulk and slab 2D-crystals are subsets of this class: bulk crystals are much 
larger in their extruded dimension than a lattice spacing, whereas, slabs are on the order 
of or less than a lattice spacing long in their extruded dimension. These structures have 
band gaps for EM radiation propagating in the plane of periodicity. Typically, slabs have a 
well-formed band gap only for polarization parallel to the periodic plane (transverse 
magnetic polarization [TM]), and a small, incomplete band gap for polarization 
perpendicular to this plane (transverse electric polarization [TE]). On the other hand, bulk 
crystals seem to possess band gaps for both polarizations. An example of a specific 2D 
geometry is the hexagonal lattice (see Fig. 1A). 

One method for constructing a 2D lattice involves stacking elongated, hollow rods in an 
array of the desired geometry, as shown in Figs. 1A and 1B. For example, a hexagonal- 
lattice photonic crystal 20 can be constructed by stacking a plurality of rods 21 in a 
triangular array. The lattice geometry in this example is hexagonal, although many other 
geometries are possible. E, B and k label the electric field vector, the magnetic field 
vector, and the wavevector, respectively. A cross indicates a vector pointing into the 
page, whereas a point indicates a vector pointing out of the page. The E, B, and k 
groupings are meant to illustrate transverse magnetic (TM) and transverse electric (TE) 
waves propagating in the TM direction (as shown) with respect to the crystal, although the 
crystal is also birefringent along other propagation directions. It has been experimentally 
observed that such a bulk, hexagonal structure composed of hollow acrylic rods (index of 
refraction 1.61) with outer diameters of 1/2 inch and inner diameters of 3/8 inch (air-filling 
fraction [AFF] = air volume / total volume = 0.60) has a fundamental TM band gap 
centered at roughly 11 GHz and a fundamental TE band gap centered at roughly 10.5 
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GHz (for propagation in the TM direction, see Fig. 1 A). The widths of the TM and TE band 
gaps are approximately 1.5 GHz and 1 GHz, respectively. For a 20-layer, bulk, hexagonal 
structure with these specifications, the field (amplitude) transmission is roughly 1/100 at 
the bottom of the TM gap, and 1/50 in the TE gap. Microscale crystals appropriate for 
optical frequencies are typically constructed using micro-fabrication techniques. 

It is possible to provide a mathematical illustration of the existence of birefringence in a 
photonic crystal. Here is a simple model that shows birefringence in a photonic crystal 
using the Kramers-Kronig relations. 

Suppose we assume for the field transmission 

f(a>)=e B( f a ^' (1) 

where w(a>) is the complex index of refraction as a function of angular frequency go, k is the 
magnitude of the wavevector in free space {k =«> ic where c is the speed of light), and d is 
the path length through the structure. The complex index of refraction can be broken up 
into its real and imaginary parts; the real part represents ordinary refraction, and the 
imaginary part represents attenuation as the wave propagates through the structure: 

r(© ) = e -ta i"Wi w e iite ["( <i, )i ib ' ^2) 

Thus, the amplitude transmission is 

If^e-*)]" . ( 3 ) 

As part of the simple model, we assume that the amplitude transmission has a simple 
periodic square-shaped form, with unity transmission everywhere except the band gap, 
which has constant transmission of T 0 . The fundamental band gap is centered at © 0 and 
extends from(o 0 -e to © 0 +e and is thus 2e in width. It is also assumed that there are 
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higher order band gaps at ma> 0 with spans from (© 0 -s) to (© 0 +e), where m is an integer 
greater than or equal to 2 and less than some reasonable cut-off value (beyond the cut-off, 
there are no additional band gaps). For simplicity, it will be assumed that this band 
structure persists until m = 10. It will also be assumed that this simple amplitude- 
transmission spectrum is polarization dependent, with T 0 , cd 0 , and e different for TM and TE 
polarizations. For the purposes of the model, the approximate specifications described 
above for the center frequencies and widths of the fundamental TM and TE band gaps of a 
20-layer, hexagonal-lattice crystal will be used. It will also be assumed that t 0 »i/ioo for 
TM waves and t 0 *i/50 for TE waves. The crystal thickness d of a 20-layer, hexagonal- 
lattice crystal with component rod outer diameter of 1/2 inch is approximately 22 cm. 

The Kramers-Kronig relations can be used to directly calculate Re[n(©)]-i given im[n(co)], 
and vice versa. In this case, we use the assumed amplitude transmission characteristics 
to determine im[n((o)] from Equation (3), then use one of the Kramers-Kronig relations to 
calculate Re[«(a>)]-i. Performing this calculation as described for TE and TM waves, a 
clear birefringence (i.e., a difference between the effective TM and TE indices of 
refraction) is found with Re[n(o>)] larger for TM waves in the transparent region below the 
fundamental band gap, as shown in the plot of Fig. 7. Similar results can be obtained in 
the transparent regions between the band gaps as well. 

This simple mathematical model is meant to demonstrate the basic principle that a 
polarization-dependent band structure is equivalent to birefringence in the refractive 
properties of the crystal. It is a simplification in that the band gaps are assumed to have 
simple square shapes and higher frequency gaps have the same widths and depths, and 
occur only at integer multiples of the fundamental. In actuality, these parameters, as well 
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as the shapes of the band gaps depend on the polarization as well as the order of the gap. 
Thus, in general, the actual birefringence can be substantially different from the results of 
this simple model. In fact, it will be shown from experimental results that the TE index of 
refraction is actually larger than that of the TM polarization for this crystal. Because of its 
simplicity, this model is not meant to provide any quantitative results, but to illustrate that 
birefringence can be inferred from a polarization-dependent transmission spectrum. 
Photonic Crystal Waveplates 

The model described above shows a relatively large birefringence (i.e., on the order of 
a few percent) in a transparent spectral region of the crystal. Since the actual 
birefringence of a 2D photonic crystal should typically be larger than this, photonic crystals 
for the first time may be utilized as compact waveplates (e.g., half- and quarter- 
waveplates) in the transparent regions between their band gaps as shown in Figs. 2A and 
2B. 

The photonic-crystal quarter-waveplate (QWP), shown generally as 40 in Fig. 2A, is a 
birefringent photonic crystal whose thickness is such that TM and TE polarized waves 
experience a relative phase shift of n/2 or 90°. A QWP will, in general, create an elliptical 
(or circular) polarization 42 from an incident linear polarization 41 and vice versa. If the 
incident linear polarization is oriented at 45° to the optic axis of the crystal, a circularly 
polarized wave is created. 

The novel photonic crystal half-waveplate (HWP), shown generally as 60 in Fig. 2B, is 
a birefringent photonic crystal whose thickness is such that TM and TE polarized waves 
experience a relative phase shift of n or 180° upon transmission. A HWP such as 60 will 
rotate the polarization of an incident linearly polarized wave 61 by -29 (if 9 is the angle 
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between the incident polarization and the optic axis of the crystal) to a transmitted wave 62 
rotated as shown in Fig. 2B. Especially at optical frequencies, these devices are 
particularly novel. Since a conventional optical waveplate usually has a birefringence 
orders of magnitude smaller than that of a photonic crystal, the thickness of a conventional 
waveplate must typically be orders of magnitude thicker than a comparable photonic 
crystal waveplate. Furthermore, since the band structure of a photonic crystal is strongly 
dependent on its geometry and material composition, its birefringent properties can be 
tailored to the particular situation and circumstances of the application. A photonic crystal 
waveplate would be used just as a conventional waveplate, with the optic axis of a 
photonic crystal waveplate is defined by the anisotropy of the structure. For example, the 
optic axis A-A of the 2D, bulk, hexagonal-lattice crystal described above, and as shown in 
Fig. 1B, is parallel to the longitudinal axes of the component rods. 

A photonic crystal waveplate of the present invention can be used to manipulate the 
polarization of transmitted light. However, photonic crystals are unique in that they can 
also be used to manipulate the polarization of reflected light. In particular, light that is not 
transmitted by a photonic crystal (i.e., blocked by the band gap) is reflected. Thus, the 
reflection spectrum of a photonic crystal must be polarization-dependent just as the 
transmission spectrum. Just as in the transmission case, this polarization-dependent 
reflection spectrum also equates with birefringence in the phases imparted to reflected 
waves. Thus, a photonic crystal can be used for the first time as a reflection waveplate in 
regions where the band gaps for the independent polarizations overlap (i.e., spectral 
regions in which both polarizations are completely reflected). In particular, it is possible to 
design photonic crystal reflection waveplates for which the reflected TM and TE 
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polarizations are completely reflected and develop a relative phase of 1/4 or 1/2 of a wave 
upon reflection. Reflection waveplates based only on dielectric materials are novel 
devices, and have no direct conventional analog. 
Photonic Crystal Polarizer 

Another device useful for polarization control is the novel linear polarizer of the present 
invention, as shown generally as 80 in Fig. 2C. A photonic crystal can be used for the first 
time as a polarizer by exploiting the separation between the band gaps for the 
independent polarizations. In particular, any crystal, which has a spectral region in which 
only one of the independent polarizations is entirely transmitted, can be used as a 
polarizer. Thus, a good photonic-crystal polarizer should have well separated band gaps. 
In the region of non-overlap, one polarization must be almost entirely transmitted, while 
the other is almost completely reflected. A 2D photonic crystal polarizer such as 80 would 
be used in the same manner as a conventional polarizer. Photonic crystal polarizers will 
be useful in that relatively thin crystals should produce high contrast in the transmission of 
different polarizations. In addition, their bandwidths and other properties can be 
manipulated by tailoring the geometry and material composition of the crystal. 
Furthermore, they do not rely on absorption, and thus, they are not heated by incident 
radiation. For example, some conventional optical polarizers absorb the power they do 
not transmit, and are therefore heated and can be easily destroyed by lasers of even 
modest power. 

Since the photonic crystal polarizer 80 interacting with incident beam 81 reflects one 
polarization shown by reflected beam 83 and transmits the other shown as transmitted 
beam 82 (Fig. 2C), it is similar to a polarizing beamsplitter. Conventional polarizing 
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beamsplitters are optical elements that produce polarized beams based on Brewster's law. 
Thus, the angle between the incident and reflected rays is limited by Brewster's angle. In 
contrast, photonic crystals do not rely on Brewster's angle and can produce polarized 
reflected and transmitted beams at any angle relative to the incident beam, depending on 
the angle between the incident beam and the surface of the crystal. The extreme example 
is a retro-reflected polarized beam, which cannot be directly produced by a conventional 
polarizing beamsplitter. 

Both slab and bulk 2D photonic crystals can be used as waveplates and polarizers. 
Although these devices have been described here so far in the context of bulk crystals, 
they can be realized with slabs as well as shown in Fig. 3. An example of a 2D photonic 
slab (side view) is shown as 100. Slab crystals have significantly different band structures 
than bulk crystals, but slab crystals should also be useful for polarization control because 
their transmission band structures are highly polarization-dependent. This design should 
be useful for planar photonic crystal circuits. As previously mentioned, slab crystals show 
only a small band gap for one of the polarizations, but a well-formed gap for the other 
polarization, and thus can be used as polarizers. Furthermore, since the degree of 
birefringence of a photonic structure is related to the degree of polarization-dependence in 
its transmission band structure, 2D slab crystals should have a very large birefringence. 
Fig. 3 shows a photonic crystal slab fabricated so that it could be used to manipulate 
polarization while still remaining flat on a planar surface (i.e., the crystal itself does not 
need to be rotated to manipulate incident horizontal and vertical polarizations); in this way, 
the optic axis can be at some angle relative to an incident polarization and the crystal can 
remain flat. The compact polarization control afforded by photonic crystals will be useful in 
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photonic crystal circuits. Photonic crystal circuits are compact optical devices based on 
photonic crystal technology. They require compact, integrable polarization control for their 
operation. 

Experimental evidence for the birefringence in photonic crystals and some potential 
applications are also shown and described herein. The data presented was collected 
using microwave radiation. Microwave transmission through a hexagonal-lattice crystal, 
as shown in Figs. 1A and 1B (see previously described specifications), was measured with 
a Hewlett-Packard vector network analyzer (VNA). Radiation was coupled to and from 
free space with microwave gain-horn antennae. The crystal was placed in a shielded box 
with an aperture (measuring 14 x 17 cm 2 ) smaller than the cross-section of the crystal, but 
large enough to minimize diffraction effects. The shielded box was placed in the far-field 
of the transmitter horn, and the receiver horn was positioned directly behind the crystal 
inside the box. The aperture and positioning were used so that the microwaves incident 
on the crystal would be effectively plane waves. The horn antennae are sensitive to only a 
single linear polarization; thus, it was possible to measure the phase delay as a function of 
frequency for microwaves propagating in the TM direction relative to the crystal, with both 
TM and TE linear polarizations. Using the measured phase delay, an effective index of 
refraction could be calculated for each polarization. It was found that there is a difference 
on the order of 10% between the indices of refraction for the two polarizations. Since the 
birefringence is a result of interference effects in the asymmetric crystal, the indices of 
refraction can depend on whether the number of layers in the crystal is even or odd. 
Representative experimental results for the TM and TE indices of refraction are shown in 
Figs. 4(a) and 4(b) for the transparent regions above and below the fundamental band gap 
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of 8- and 16-layer hexagonal-lattice crystals, respectively. Using the measured 
birefringence, it is possible to compute the frequencies at which TE and TM waves 
propagating through the same length of crystal would accumulate a net phase difference 
of 1/4 or 1/2 of 2n. At these frequencies, the crystal acts as a quarter- or a half- 

waveplate, respectively. Specifically, if the total relative phase shift is |(2m+i) for an 

integer m greater or equal to 0, the crystal acts as a quarter-waveplate of (m+i)' A order. 
On the other hand, a relative phase shift of n(2m+i) is produced by a half-waveplate of 
{m+if order. Thus, for a quarter-waveplate of (m+if order at frequency w, the 
birefringence must be 



In what follows, the right-hand sides of Equations (4) and (5) will be defined as the quarter- 
wave condition (QWC) and half-wave condition (HWC), respectively, of (m+if order. 

In Figs. 4C and 4D, the differences between the indices of refraction for 8- and 16-layer 
crystals along with calculated QWCs and HWCs of the indicated order are shown. The 
intersection of a quarter- or half-wave condition curve with the experimentally measured 
index difference indicates that the crystal acts as the corresponding waveplate at that 
frequency. 

Using these experimental predictions, 8-layer quarter- and half-waveplates were 
tested. By rotating the crystal (about an axis parallel to the TM direction) so that the 




(4) 



Similarly, for a half-waveplate, we have 




(5) 
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incident polarization makes an angle of 45° with the optic axis (here, the optic axis is 
parallel to the component rods and the incident wave propagates completely in the 
periodic plane), a circularly polarized wave should be produced if the incident wave has 
the quarter-waveplate frequency. At this frequency, the receiver horn should detect the 
same power transmission as it is rotated to any angle about the TM direction (the receiver 
horn detects the linear projection of the circular polarization). The results of this 
experiment with an 8-layer crystal are shown in Fig. 5A. Notice that all the traces cross at 
the same frequency predicted by the data in Fig. 4C. 

In another set of experiments, both horns were rotated simultaneously (about an axis 
parallel to the TM direction) in the aligned configuration, while the transmission at a fixed 
frequency was measured as a function of the receiver horn angle (horns parallel). This 
procedure was also repeated with the transmitter-horn antenna maintained at a constant 
90° relative to the receiver-horn antenna as both were rotated simultaneously (horn 
antennae crossed (i.e., perpendicular)). For a quarter-waveplate, the expected angular 
dependencies are 

t„ = Vcos 4 9+sin 4 e (6) 

and 

'i=^tM29) (7) 

for the normalized transmitted field amplitudes in the parallel and perpendicular 
configurations, respectively, where 9 is the receiver horn angle (defined to be zero when 
oriented in the direction of a TM polarization). On the other hand, for a half-waveplate, the 
expected angular dependencies are 



18 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 



7},= cos 2 (29) (8) 

and 

71= sin 2 (20) (9) 

for the normalized power transmission in the parallel and perpendicular configurations, 
respectively. 

For both quarter- and half-waveplates, the measured results at the waveplate 
frequencies agree very well with the expected angular dependencies for the corresponding 
waveplates (see Figs. 5B and 5C). For the quarter-waveplate, the data was obtained by 
measuring the normalized field transmission at approximately 17.4 GHz; for the half- 
waveplate, the normalized intensity transmission was averaged over 15-15.5 GHz. These 
results clearly show, for the first time, the successful design and operation of compact, 
photonic-crystal waveplates. 

Experimental verification of the use of a photonic crystal as a polarizer or polarizing 
beamsplitter is also described herein. Using the same experimental setup described 
above, the transmission through a 20-layer, hexagonal-lattice crystal was measured with 
the horn antennae parallel and perpendicular (as a function of receiver horn angle) at a 
fixed frequency within the band gap of one polarization, but in the transparent region of the 
other. The results of this trial are shown in Fig. 6A for the fundamental band gap. This 
experiment was also performed utilizing the next order band gap of the structure (see Fig. 
6B). The main portions of these figures show the field transmission profiles of the crystal 
and the insets show the results of the horn rotation operations at roughly 1 1 .4 GHz (a) and 
19.2 GHz (b). In both cases, the experimental results agree quite well with the expected 
angular dependencies for these operations. 
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Higher order band gaps occur at multiples of the fundamental band gap frequency. In 
addition, other band gaps exist at frequencies higher than that of the fundamental. Since 
band gaps at higher frequencies can be wider and/or deeper than the fundamental band 
gap, it is possible for the operation of a photonic crystal polarizer to be more effective at 
higher frequency; thus, at fixed frequencies, it may be advantageous to scale down the 
crystal in size to utilize higher order band gaps. Since the photonic crystal polarizer 
reflects the light that is not transmitted, the reflected polarization can be utilized, as well. 
These results clearly show directly, for the first time, the successful operation of a photonic 
crystal polarizer/polarizing beamsplitter. 
Photonic Crystal Reflection Waveplate 

Another useful property for polarization control is the polarization-dependent phase 
shift acquired upon reflection from a photonic crystal. A polarization-dependent phase 
shift may be acquired upon reflection from (or transmission through) a photonic crystal 
even at normal incidence if the crystal is periodic in more than one dimension. For 
example, in a crystal with two-dimensional (2D) periodicity, the Fresnel coefficients for the 
two basis polarizations may be different at each material interface, and the cumulative 
effect of many curved interfaces leads to a polarization-dependent band structure. These 
polarization-dependent effects can also occur for light obliquely incident on a one- 
dimensional (1D) crystal. 

A photonic-crystal-based reflection waveplate can be constructed completely from 
dielectric materials (although other components are not precluded) and it does not require 
external mirrors or additional optical elements. The reflective and polarization-dependent 
properties of a photonic crystal can be tailored by changing the lattice geometry, material 
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composition, utilizing higher order band gaps, and changing the incident angle of the 
radiation; thus, photonic crystal waveplates (both reflection and transmission) can be 
adapted for use in many possible polarization-control applications. Furthermore, they can 
be used at any wavelength scale simply by applying the appropriate scaling. An 
experimental proof-of-principle demonstration is presented here. 

For a practical reflection waveplate, this polarization-dependent phase delay is most 
interesting in spectral regions of high reflectivity, i.e., where both polarizations are almost 
completely reflected. The physics behind a polarization-dependent phase shift is again 
based on the cumulative effect of the complex Fresnel coefficients at each curved material 
interface; however, in the band gap portion of the spectrum, material interfaces deeper 
into the crystal contribute increasingly less phase shift to the reflected wave because the 
field is exponentially attenuated within the crystal. 

Unfortunately, knowledge of the phase delay of a transmitted wave in transparent 
spectral regions cannot be used to determine the phase of the reflected component in the 
band gap. Thus, the phase delay experienced by reflected waves in the band gap 
remains an open question, despite phase measurements of photonic crystal transmission. 
Here it is demonstrated that the relative phase between the two basis polarizations 
acquired upon reflection from a 2D photonic crystal can be significant in spectral regions of 
very high reflectivity. This property is applied to the construction of a microwave-scale, 
reflection quarter-waveplate based on the 2D, square-lattice geometry and its operation is 
demonstrated. Since the properties of photonic crystals are based on the physics of 
Maxwell's equations with boundary conditions, the results presented here are scale- 
invariant and apply across the electromagnetic spectrum. 
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The crystals were constructed, as shown schematically in Fig. 8, by stacking acrylic 
(index of refraction 1.61) pipes 121 in an array of the desired geometry (e.g., square 
lattice) and cementing them together with an acrylic solvent. The crystal used for this 
experimental demonstration was a bulk (i.e., the length of the extruded dimension was 
much longer than a lattice constant), 30-layer, square-lattice structure, with a lattice 
constant of 1/2 inch, and an air-filling fraction (AFF) of 0.66. The square-lattice 
arrangement is shown in Fig. 8, forming photonic crystal 120. 

Transmission measurements of the fundamental band gaps of this crystal were made 
with a Hewlett-Packard vector network analyzer (VNA) using the methods previously 
described. Radiation was coupled to and from free space with microwave gain-hom 
antennae. The crystal was placed in a shielded box with an aperture smaller than the 
cross-section of the crystal. The shielded box was placed in the far-field of the transmitter 
horn, and the receiver horn was positioned directly behind the crystal inside the box. The 
aperture and positioning were used to eliminate edge effects from the crystal, and to 
ensure that the microwaves incident on the crystal would be effectively plane waves. The 
horn antennae are sensitive to only a single linear polarization; thus, it was possible to 
measure the phase delay as a function of frequency for microwaves propagating in the 
periodic plane of the crystal, with both TM and TE linear polarizations (we label the 
polarizations with respect to the periodic plane: the transverse magnetic (TM) and 
transverse electric (TE) polarizations are parallel and perpendicular to the periodic plane, 
respectively (see Fig. 1). Since the crystal is non-absorbing, a transmission measurement 
suffices to show the power reflectivity of the crystal. The experimentally measured power 
reflectivity of the crystal as a function of frequency in the vicinity of the fundamental band 
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gaps is shown in Fig. 9. From these transmission measurements, we find that the 30-layer 
crystal exhibits greater than or equal to 95% power reflectivity for both TM and TE 
polarizations over roughly a 100 MHz overlap range centered near 9.6 GHz. 

Reflection measurements were also made using the experimental setup presented in 
Fig. 8. Microwaves were coupled between the vector network analyzer (VNA) 130 to and 
from free space with polarization sensitive horn antennae. Transmitter 150 and receiver 
160 horns were mounted together on a rotation stage inside a microwave-shielded box 
140, internally lined with microwave absorbing material (Eccosorb). The crystal was 
placed in the far-field of these horns, behind an aperture 141 smaller than the cross- 
section of the crystal 120. Waves emitted from the transmitter horn 150 were incident on 
the crystal and retro-reflected back to the receiver horn 160 juxtaposed with the 
transmitter. The rotation stage allowed the horns to be rotated simultaneously in the plane 
parallel to the incident face of the crystal. The horns could be mounted with their 
antennae either parallel or crossed to emit and detect either parallel or orthogonal 
polarizations, respectively. The reflection measurements show that the phase difference 
between TM and TE polarizations upon reflection from this crystal is approximately 90° 
over a continuous span of roughly 100 MHz just below 9.6 GHz (see Fig. 10). These data 
suggest that this photonic crystal can be used as a reflective quarter-waveplate with high 
reflectivity in the spectral region where its band gaps overlap. 

This hypothesis was confirmed with direct reflection measurements showing the 
behavior of a quarter waveplate. Fig. 11 shows the experimentally measured field 
reflection from the 30-layer square lattice crystal, at a fixed frequency just below 9.6 GHz, 
as a function of transmitter horn angle with the horns parallel (open circles) and crossed 
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(filled squares). The data were obtained from raw measurements as follows: the results of 
multiple measurements at the same angle were averaged to wash out alignment 
uncertainties in the crystal and horns (system was realigned from scratch for each 
measurement); an offset was added to the points belonging to the horns parallel data set, 
and a normalization factor was applied to the horns crossed data set. The offset and 
normalization factor were used to bring the data sets to their expected maxima, and were 
needed because of the inherent difficulties in calibrating free-space microwave reflection 
measurements to an absolute amplitude scale. The resulting data sets agree very well 
with the expected theoretical curves also shown in Fig. 1 1 . 

In addition to the square lattice, some other bulk structures were also constructed and 
tested. However, for these other structures, the maximum phase difference between 
reflected TM and TE polarizations was observed to be less than 90° over nearly all of the 
overlap region of the fundamental band gaps. These results stress the importance of the 
lattice properties in the design of photonic crystal devices for polarization control. 

In the experiments discussed here, normally incident radiation was considered, 
although other angles of incidence are also permissible. For non-normal incidence, 
however, the waveplate will have a modified frequency range of operation. If the incident 
beam is obliquely incident on the crystal, the angle of the reflected beam will be 
determined by the angle of incidence. Thus, in principle, the angle between the incident 
and the reflected beam can be anything from 0° to 180°. In addition, shifts in incidence 
angle are known to shift band gap position, and should modify phase-dependent 
properties, providing another degree of freedom to tailor the properties of the waveplate. 
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These experiments utilized the fundamental band gaps of the photonic crystal. Higher 
order band gaps occur at integer multiples of the fundamental band gap frequency. Since 
higher frequency band gaps can be wider and deeper, reflection waveplates operating in 
these higher order gaps can potentially have greater reflectivity and bandwidth. 

In conclusion, this is a proof-of-principle experimental demonstration of an all-dielectric, 
reflection waveplate based on a photonic crystal. It is possible to tailor the properties and 
improve the characteristics of these waveplates by modifying lattice geometry, material 
composition, filling fractions, angle of incidence, and utilizing higher-order band gaps. If 
the incident radiation is normally incident on the face of the crystal, it is retro-reflected; 
however, if the incident beam is obliquely incident on the crystal, the angle of the reflected 
beam is determined by the angle of incidence. Thus, in principle, the angle between the 
incident and the reflected beam can be anything from 0° to 180°. Finally, as seen in Fig. 
10, the TM-TE phase difference can well exceed 90°, indicating that reflection quarter- 
waveplates are not a limiting case; with certain crystal configurations, reflection half- 
waveplates should also be possible. Dielectric, reflection waveplates may have important 
practical applications, particularly in optics, where compact, designable polarization control 
is desirable for integrated, photonic devices. 

A Wide-Band ( Achromatic) Photonic Crystal Transmission Waveplate 

A wide-band (achromatic) transmission waveplate can be constructed if the photonic 
crystal is designed such that the desired operating wavelength is larger than that of the 
fundamental band gaps. In this spectral region, the crystal is transparent, and the 
birefringence does not change significantly with frequency. These types of waveplates 
can be much more compact than conventional wide-band polarization control devices, and 
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their specific dispersive properties and bandwidth can be tailored by modifying the crystal 
geometry and material composition. The principles behind polarization control with 
photonic crystals can be applied at any wavelength simply by applying the appropriate 
scaling. These devices have many potential applications; for example, they should 
provide highly integrable polarization control in optical photonic circuits. Polarization- 
dependent effects in the achromatic photonic crystal transmission waveplate can also be 
understood in terms of the cumulative effect of polarization-dependent Fresnel coefficients 
at each material interface. However, in spectral regions well below the fundamental band 
gap, the effective indices of refraction are non-dispersive (i.e., not frequency-dependent). 
This property can be qualitatively understood by considering the Kramers-Kronig relations 
applied to the transmission model previously presented (Fig. 7). From this simple model, 
there is a clear birefringence that becomes relatively non-dispersive at low frequency. 

Experimental evidence for the birefringent behavior characteristic of a wide-band 
waveplate is also presented here. The crystal used was a bulk, hexagonal-lattice structure 
composed of stacked acrylic pipes (lattice constant 1/4 inch, AFF = 0.32). The data 
presented were collected using the experimental transmission setup described above. 
Fig. 12 shows representative results for the experimentally measured birefringence 
between TM and TE waves transmitted through a 20-layer crystal. Experiments show that 
crystals with different (even) numbers of layers exhibit similar birefringences, but different 
phase retardances owing to the different path lengths. As expected, the birefringence is 
essentially constant over measured interval. This fact is dramatic considering that the 
interval covers a large frequency range: the ratio of the maximum to minimum frequency is 
more than 2.5. Furthermore, the size of this interval was only limited by the frequency 
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range of the measurement apparatus. Also shown in Fig. 12, is the quarter-waveplate 
condition, illustrating that this crystal acts as a quarter-waveplate with a scaled useful 
bandwidth comparable to standard, commercial, zero-order waveplates for optical 
wavelengths. 

Although the measured birefringence is nearly linear, it evident from the data that it has 
local frequency variations or wiggles that become more pronounced at wavelengths that 
are larger compared with the lattice spacing. At first glance, these variations may seem to 
be measurement noise; however, using a simple 1 D-photonic-crystal model, it can be 
shown that they are in fact an inherent property of photonic crystals. Although 1D 
photonic crystals are not birefringent at normal incidence, it is still useful to study their 
dispersive properties from their effective or geometrical index of refraction (calculated from 
the transmitted phase). This model displays nontrivial features, many of which can be 
applied, at least qualitatively, to the design of birefringent, 2D structures. 

Fig. 13 shows the calculated effective indices of refraction below the fundamental band 
gap of three different 1D, acrylic-air stacks of 20, 40, 60, 80, and 100 layers (plotted as a 
function of frequency, normalized to the central frequency of the fundamental band gap of 
each structure). One stack consisted of quarter-wave layers, and thus had an AFF of 
roughly 0.62. By reducing the thickness of the air layer and increasing that of the acrylic 
layer by the same amount, or vice versa, stacks of two other AFFs (0.32 and 0.92) were 
created. For simplicity, the raw materials were assumed to have zero intrinsic dispersion. 
These calculations show several interesting features. First of all, we see that the indices 
of refraction are flat overall, but begin to increase as the band gap is approached; 
however, the slope of the increase near the gap apparently depends on the AFF, with the 
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quarter-wave stack showing a greater slope than the other two. Since the quarter-wave 
stack has the widest, deepest band gap of the three structures, the band gap width and 
depth are evidently important features for the dispersion at higher frequencies. Secondly, 
it is clear that this model exhibits local frequency variations in the index of refraction, which 
were also seen in the experimental data of Fig. 12. As seen in the experimental data, the 
amplitude of these variations becomes more pronounced at lower frequency, but the 
amplitude is reduced and the rapidity increased as more layers are added. In addition, the 
variations show a clear correlation with AFF, becoming less significant as the AFF is 
increased. 

These calculations suggest that photonic crystal waveplates are highly engineerable 
devices. Furthermore, the local variations in the indices of refraction can be exploited to 
create a waveplate in which the relative phase retardance is virtually frequency- 
independent over some frequency range. In particular, the variations exhibit anomalous 
dispersion (i.e., over some frequencies, the index of refraction decreases with increasing 
frequency); if the birefringence is designed to decrease as roughly as (or nearly so) 
over some interval, the phase retardance will be constant over that range. 

It has been shown here that it is possible to control the polarization of EM waves 
propagating in the plane of periodicity of a two-dimensional (2D) photonic crystal (i.e., 
photonic crystals periodic in two dimensions). The principles of polarization control with 
photonic crystals described here are also valid for other classes of anisotropic photonic 
crystals; for example, three-dimensional (3D) photonic crystals (i.e., photonic crystals 
periodic in three dimensions) can also have the polarization-dependent boundary 
conditions necessary for polarization control. 
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Nonlinear Optical Frequency Generation with a Photonic Crystal 

Nonlinear optical crystals can be used to generate new frequency components from 
input light beams. In such interactions, phase-matching conditions must be satisfied for 
efficient conversion to take place. This is usually accomplished by angle tuning 
birefringent crystals to create the correct indices of refraction for ordinary and 
extraordinary rays involved in the nonlinear mixing process in order to phase-match. 
However, the drawback of this method is the walk-off or divergence of the ordinary and 
extraordinary waves as they propagate through the crystal, an effect that limits the 
efficiency of the nonlinear mixing process. If instead, a birefringent photonic crystal 
composed of material with the correct optical nonlinearity is used as the nonlinear optical 
medium, the birefringence can be designed to satisfy the phase-matching conditions at the 
desired frequencies by tailoring the lattice geometry and size, and the material proportions 
and composition. Then, phase matching can be achieved without the need for angle 
tuning, increasing the efficiency of the nonlinear mixing process. On the other hand, it 
may also be possible to achieve phase matching for waves of a single polarization by 
combining the effects of material and geometrical dispersion. 

Birefringent photonic crystals composed of material with optical nonlinearity can be 
designed to increase the efficiency of nonlinear mixing processes. As a specific example, 
consider x (2) processes, which can occur in non-centrosymmetric crystals. In this type of 
process, it is possible to produce sum-frequency, second-harmonic, and difference- 
frequency (parametric amplification) generation. 

In order for efficient conversion to take place, conservation of momentum between the 
waves involved in the nonlinear mixing process must be satisfied. In particular, the 
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nonlinear mixing process is efficient if the wavevector mismatch M is zero. Stated 
another way, the waves must have the correct relative phasing in order for energy to be 
efficiently extracted from the input beam. This condition is known as phase matching. For 
example, in the nonlinear process of second-harmonic generation, an input beam with 
wavevector ki and frequency <d, can be converted to an output beam with wavevector 
and frequency ©2=2(0,; in this situation, the phase matching condition is Ak = 2k x - k 2 = 0 
for collinear beams. This condition is equivalent to the requirement n(<o,)=n(2a>,) where a>, 
is the frequency of the fundamental harmonic and «(co ) is the frequency-dependent index 
of refraction of the optical medium. For collinear m th harmonic generation (for positive 
integer m >1), wherein an input beam at frequency a>, is converted to an ouput beam at 
frequency mco, in a medium with % {n) nonlinearity, phase matching requires 
Ak - mk x -k 2 = 0 , which is equivalent to «(©,)= n(wto, ). 

For many nonlinear wave mixing processes, this phase matching is normally achieved 
in the prior art by angle tuning a birefringent nonlinear medium. For example, in second 
harmonic generation, the angle of the medium is adjusted such that an ordinary polarized 
wave of the fundamental harmonic and an extraordinary wave of the second harmonic (or 
vice versa) satisfy this condition. This is possible because the index of refraction for 
extraordinary polarized waves is a function of the angle (9) of propagation relative to the 
optic axis of the medium. However, the Poynting vectors (i.e., the direction of energy 
transport) of the ordinary and extraordinary rays will not be parallel whenever the 
extraordinary ray propagates at an angle that is not 6 = 0° or 9 = 90° relative to the optic 
axis. As a result, the spatial overlap of the ordinary and extraordinary polarized beams will 
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be reduced, decreasing the efficiency of the nonlinear process, whenever 0 is not 0° or 
90°. A birefringent photonic crystal composed of material with the x (2) optical nonlinearity 
can be designed so that the ordinary and extraordinary waves satisfy this condition at 0 = 
0° or 0 = 90°, eliminating the walk-off of the ordinary and extraordinary waves. 

Furthermore, many materials with x (2) nonlinearity are not intrinsically birefringent, 
preventing them from being used in any phase-matching schemes that require 
birefringence. However, a birefringent, photonic crystal structure can be constructed using 
these materials; as described above, a photonic crystal can be birefringent due to its 
geometrical properties, even though the component materials may have no intrinsic 
birefringence. Thus, such a nonlinear material can be used in birefringent phase-matching 
schemes if it is a component in a birefringent photonic crystal, even if the material has no 
intrinsic birefringence. 

In some cases, another established process known as temperature tuning can be used 
to achieve phase matching for ordinary and extraordinary rays without deviating from 0 = 
0° or 0 = 90°. This scheme exploits the temperature dependence of the birefringence of 
certain crystals. However, use of a birefringent photonic crystal provides a novel 
alternative, which may be especially useful for nonlinear media whose birefringence is not 
significantly temperature dependent. 

On the other hand, photonic crystals can also offer another method of phase matching 
which does not rely on any birefringence at all. Usually, normal dispersion prevents the 
phase matching condition from being satisfied for two waves of the same polarization 
since they are at very different frequencies. However, the structural and compositional 
properties of photonic crystals produce an effective geometrical index of refraction. The 
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dispersive properties of this geometrical index can be tailored by tuning to different 
frequencies above or below the band gaps of the media. Thus, for example, it is possible 
to achieve phase matching for second-harmonic generation by designing a structure for 
which material dispersion and geometrical dispersion cancel or combine in the right way to 
make «(o> , ) = /i(2o> , ) . 

Second-harmonic generation provides an example of the utility of a photonic crystal for 
increasing the efficiency of nonlinear optical mixing process; however, it is not a limiting 
case - the birefringent and/or dispersive properties of photonic crystals can be tailored to 
suit other nonlinear optical mixing processes as well. 

Fig. 14 shows a schematic of a nonlinear photonic crystal used for second harmonic 
generation. 
Delay Lines 

An experimental and theoretical study of a simple, passive system consisting of a 
birefringent, two-dimensional photonic crystal and a polarizer in series is presented. It is 
shown that fast and slow group velocities can arise even though no incident radiation is 
absorbed or reflected. This type of system is an example of a new, broad class of 
systems that can be used as delay lines for electromagnetic waves and electronic signals; 
the central characteristic of these systems is a coupling between two or more modes. A 
vector formulation of the Kramers-Kronig dispersion relations facilitates an understanding 
of these counter-intuitive effects. 

Superluminal group velocities have been observed in a number of different physical 
systems. These include passive absorptive, passive reflective, and active transparent 
media. Here the first discussion and experimental observation of superluminal effects in a 



32 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 



passive system with neither absorption nor reflection is presented. The effects arise 
because of a transfer of energy or interference between two modes of the electromagnetic 
field, in this case two different polarizations of light. These results can be interpreted using 
a new vector formulation of the Kramers-Kronig relations. Using this formalism, it is clear 
that these arguments can be applied to any systems (both passive and active) in which 
energy can be transferred between two or more modes. 

It is widely believed that the Kramers-Kronig (K-K) relations require that passive 
systems be either absorptive or reflective in order to exhibit superluminal effects. 
However, this is not the case; while absorptive and reflective systems are required to have 
spectral regions of anomalous dispersion, the converse is not necessarily true. In fact, the 
exchange of energy between modes is a sufficient condition for anomalous dispersion in 
any system. These effects are consistent with causality. 

These effects can be observed in a broad range of systems. As a specific example, 
consider an experimental system consisting of a slab of highly birefringent, 2D photonic 
crystal and a linear polarizer, placed in series. The role of the photonic crystal here is its 
large, frequency-dependent birefringence in transparent spectral regions; however, a 
conventional birefringent material could potentially be used instead. The results presented 
here were obtained using a photonic crystal with fundamental and second-order photonic 
band gaps in the regions of 10 and 20 GHz, and displaying strong birefringence with very 
high transmission in the frequency range between the two gaps. The crystal itself was a 
bulk, 18-layer, hexagonal-lattice array of hollow acrylic rods (outer diameter 1/2") with an 
air-filling fraction (AFF) of 0.60. 
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The transmission and dispersive properties of this system between the two band gaps 
were studied using a Hewlett-Packard vector network analyzer (VNA). Microwaves were 
coupled to and from free space with polarization-sensitive horn antennae. The photonic 
crystal was placed in the far-field of the transmitter horn to ensure that planar wavefronts 
of a well-defined polarization were incident on it. The receiver horn was positioned 
immediately behind the photonic crystal on a direct line of sight with the transmitter horn. 
In addition, the crystal and receiver horn were placed inside a microwave-shielded box 
with an open square aperture, whose size was chosen to minimize diffraction effects while 
eliminating signal leakage around the crystal. 

In order to control the polarizations of the incident and detected fields relative to the 
fast axis of the crystal, the transmitter and receiver horns were mounted on precision 
rotation stages. The angle of the incident polarization 9 was held fixed at 45° relative to 
the fast axis of the crystal, while the angle of the receiver horn p was allowed to vary. The 
coordinates are defined such that the slab is oriented with its fast axis parallel to the 
vertical direction, and label the incident polarizations TM (transverse magnetic) and TE 
(transverse electric) for polarizations parallel and perpendicular to the fast axis, 
respectively. The experimental setup is displayed in Fig. 15. The transmission amplitude 
and phase measurements with the receiver horn set at 40° and 50° are shown in Figs. 16A 
and 16B, respectively; calculated (dotted lines) and measured (solid lines) transmission 
(left axis) and phase (right axis) for the detected polarization inclined at (a) 40° and (b) 50° 
are shown. Since the spectral region shown is far from the band gaps, the observed 
transmission dip is not due to any band gap effect. It is caused by the fact that the 
photonic crystal rotates the polarization of the light by adding a different frequency- 
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dependent phase to each polarization component. For p<45° there is clear anomalous 
dispersion in the vicinity of 16.5-17 GHz (the half-waveplate frequency for the photonic 
crystal) while for p > 45° the dispersion is normal. A remarkable feature of this data is that 
while the transmission is identical for both p=40° and p=50°, aside from experimental 
errors, the phase properties are quite different. These features can be explained using a 
physical model that is developed below. 

A linear system that is invariant under time translation and is described by a scalar 
response (Green's) function g (t) produces a time-dependent response b(t) to an input a(t) 

given by the convolution b{t)= £g(t-T>(t)rfT . This expression is the starting point for the 

usual derivation of the Kramers-Kronig relations. The present system, however, accepts a 
time-dependent vector input and produces a vector output. The scalar convolution must 
therefore be replaced with the expression 

*■(<)= f Sy{t-^^)dx = P gyit'^jit-ttyS (10) 

*-00 J- CO 

where b,{t) and a,{t) are the /'* components of the time-dependent output and input 

vectors b(*) and *{t), and g „(t) is the (ijf component of the Green's function matrix which 

describes the system. The standard summation convention has been employed over 
repeated indices. In the frequency domain, Equation (10) takes the simple form 

BfahGy^Ajfa) (11) 

where £,(<»), ^(o), and g ? (©) are the Fourier transforms of b,(t), a ( {t), and gg (t), 
respectively (assuming these transforms exist). 
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Causality dictates that the output must vanish for times before the input has 
propagated through the system. Thus, we require £ g^^t-x^x =o where T must be 

greater than or equal to the relativistic propagation time. Since this relation must be 
satisfied for any choice of input, the components of the Green's function matrix vanish 
individually for all time prior to T. This directly implies that the components G i , (a>)r' m7 ' have 
analytic continuations for complex frequencies $ where im(<o)>o (i.e., the upper half-plane 
in complex frequency space). Therefore, it is possible to proceed in the usual way to show 
that the real and imaginary parts of each component g^) satisfy the Kramers-Kronig 

dispersion relations assuming the components G l7 (a>) are square integrable: 



where P denotes Cauchy's principal value. 

The signals in our experiment can be expressed as A(©)=^( ( o)e /< and B((o)=£((o> fil 
where and e B are unit vectors in real space, and 3(co) and £(o>)are the complex 
amplitudes of the incident and detected electric fields. In this special case where e A and 
e B are frequency-independent, Equation (11) implies that a scalar transfer function 
describes the relationship between the complex functions «(©) and a(<s>) 




(12) 




(13) 



4(b) 



(14) 
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where e\ and e' B are the /'* components of the unit vectors e A and e a , respectively. Since 
the elements G y (<o) satisfy the Kramers-Kronig relations, it follows that the scalar function 
#(©) must also satisfy them. However, the components of g (> ((o) interfere with each other 

in #(©). 

It is this interference which leads to the superluminal effects in the experiment 
described here. The Green's-f unction matrix for the present birefringent system is 



G(©)< 



0 



(15) 



where ^ rE (a)=n m ((o)i>d/c and § TM {a>)=n m (a)s>dic are the frequency-dependent phases 
imparted to TE and TM polarizations, and n m {<») represent the frequency- 

dependent indices of refraction for the two polarizations, d is the slab thickness, and c is 
the vacuum speed of light. Substituting these matrix elements into Equation (14), it is 
found that 

^(co) = sin(p)sin(9y*^ a> ) + cos(p)cos(ey*™ (o) . (16) 

For incident linear polarization at an angle of e =« /4 from the vertical, the magnitude and 
phase of the transfer function are, respectively, 

^(co)| = -^{l + sin(2p)cos[A<l,(«,)]} 1/2 (17) 



arg[#((o )]= arg^'*™( <,> )]+ arctan 

where A<Kffl)=<|> re (a> )-«>„,(©). 



sin[A<|>(co)] 
cos[A(|>(a))]+cot(p) 



(18) 



The transmission |#(©)| has a minimum for A<|>(co m )=(27M+i)i where co m is a half-waveplate 



frequency of the photonic crystal and m is an integer. For p=*/4, #(co m )| = o since the 
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polarization which emerges from the slab has zero projection along the unit vector e a . At this 
point, the phase arg[#((o)] is undefined since the transmission vanishes. For p =* /4±e where 
o<e «i, we find 

|//(<o m )| = e + &(e 3 ) (19) 
5arg[//((o)| _ + 1 5A<t>(co) 



3co j 2e 9© 



(20) 



where x g = 3arg|/?(© )]/&<> is the group delay of the transmitted wave. As e ->o (i.e., as we 
approach the singularity in the transfer function), the transmission goes to zero at © m and 
the group delay is unbounded. For example, if p=7i/4-e and dA^©)/^ >o it is clear 

that the group delay can become superluminal, zero, or even negative depending on the 
value of e. On the other hand, the group velocity can be arbitrarily slowed if p =n /4+t . 
Thus, if one were to measure the time of flight of an analytic pulse based on the arrival of 
its peak, any result can be obtained. However, there is no violation of causality here: the 
group velocity of the pulse has nothing to do with the signal velocity (i.e., the velocity of 
"information"), the quantity restricted by relativity. As we will demonstrate empirically, 
these superluminal group velocities are in fact required by causality through the Kramers- 
Kronig relations. 

In Figs. 16A and 16B, the results of this simple model described by Equations (17) and 
(18) are compared with the experimental results. The frequency-dependent phases $ TE {a>) 
and <t> w (a>) were calculated using an independent measurement of the indices of refraction 
of the photonic crystal. It is clear that the model correctly predicts the behavior of this 
system. Some discrepancies between the model and the experiment are observed at the 
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edges of the plotted regions due to the effects of the band structure of the photonic crystal, 
which were not included in the model. 

The interference between terms in h(g>) arises because of the two (polarization) 
modes available to light in this system. The birefringent photonic crystal allows a coupling 
between these modes, and energy can flow from one to the other. 

In the system discussed here, the propagation speed (i.e., group velocity) of a 
polarized light pulse (i.e., wavepacket) or waveform displays both fast (superluminal) and 
slow light behavior with a sharp transition between the two regimes due to the presence of 
a singularity in the transfer function of the system. This allows for total control over the 
speed of a wavepacket or waveform by adjusting the either the relative angular 
orientations of the birefringent crystal, the polarizer, and/or the incident EM wave 
polarization, and may have several important applications in optical signal processing and 
optical memory. In this particular system, the transmission goes to zero near the 
singularity; however, this is not a fundamental limitation and active systems with gain 
rather than loss can also be designed as delay lines using this principle of energy transfer 
between modes. 

The foregoing description of the invention has been presented for purposes of 
illustration and description and is not intended to be exhaustive or to limit the invention to 
the precise form disclosed. Modifications and variations are possible in light of the above 
teaching. The embodiments were chosen and described to best explain the principles of 
the invention and its practical application to thereby enable others skilled in the art to best 
use the invention in various embodiments and with various modifications suited to the 
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particular use contemplated. The scope of the invention is to be defined by the following 
claims. 
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